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1. INTRODUCTION 

How much ice is there in the Tropical 
Tropopause layer, globally? How does one 
begin to answer that question? Clouds are 
currently the largest source of uncertainty in 
climate models, and the ice water content 
(IWC) of cold cirrus clouds is needed to 
understand the total water and radiation 
budgets of the upper troposphere and lower 
stratosphere (UT/LS). The Cloud-Aerosol 
Lidar and Infrared Pathfinder Satellite 
Observation (CALIPSO) satellite, originally a 
"pathfinder" mission only expected to last for 
three years, has now been operational for 
more than eight years. Lidar data from 
CALIPSO can provide information about 
how IWC is vertically distributed in the 
UT/LS, and about inter-annual variability and 
seasonal changes in cloud ice. However, 
cloud IWC is difficult to measure accurately 
with either remote or in situ instruments 
because IWC from cold cirrus clouds is 
derived from the particle cross-sectional 
area or visible extinction coefficient. 
Assumptions must be made about the 
relationship between the area, volume and 
density of ice particles with various crystal 
habits. Recently there have been numerous 
aircraft field campaigns providing detailed 
information about cirrus ice water content 
from cloud probes. This presentation 
evaluates the assumptions made when 
creating the Cloud-Aerosol Lidar with 
Orthogonal Polarization (CALIOP) global 
IWC data set, using recently reanalyzed 
aircraft particle probe measurements of very 
cold, thin TTL cirrus from the 2006 CR-AVE. 
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2. CLOUD ICE WATER CONTENT 
PARAMETERIZATION 


Cloud ice water content is created in the 
CALIOP data set using a two-step process. 
As an elastic backscatter lidar, CALIOP 
measures attenuated backscatter, and 
extinction is retrieved using the method of 
Young and Vaughan (2009). Retrieving 
extinction for thin subvisible cirrus requires 
the assumption of extinction to backscatter 
coefficient ratio, currently 25 sr used in the 
CALIOP Version 3 algorithm. Given the 
retrieved extinctions, Version 3 uses an 
empirically based power law: 

IWC = ci<j b , a=1 19, b=1.22 (1) 

This power law is based on a set of in situ 
and remote aircraft-based measurements 
(Heymsfield et. al., 2005). 


Figure 1 illustrates these steps for an 
overpass of Hurricane Sandy, from total 532 
nm backscatter, to extinction and then IWC. 
(Note - The differing tropical (LHS) and 
extra-tropical (RHS) ice cloud morphology of 
the Tropical and extratropical sections of 
Sandy is evident in this series of plots.) 



Figure 1: Representation of the CALIOP ice 
water content parameterization process. 





3. ICE PARTICLES IN THE TTL 

Ice clouds in the Tropopause Transition 
Layer (TTL) are difficult to measure because 
they are thin and most are sub-visible. The 
TTL is the region where vertical motion 
transitions between convection and the 
large-scale Hadley circulation. Figure 2 a, b 
and c show the subtle boundary between 
the level of neutral buoyancy with maximum 
convective outflow and the TTL at 14.5 km 
(in a) or 355 K potential temperature (in b). 
These plots show data from a representative 
month (December 2009). 
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Figure 2,a-c: Occurrence of Cold Ice Clouds 
during December 2009. Shown are a) 
Potential temperature (theta) vs. mid-cloud 
altitude, b) Zonal cloud occurrence profile 


(theta y-axis), c) Integrated attenuated 
backscatter in the tropics as a function of 
altitude (theta). 

As shown in 2c, most TTL cirrus are quite 
thin. Volume depolarization greater than 
0.05 indicates that these layers are not 
aerosols. Figures 3a and 3b show the 
distribution of ice water path, depolarization 
and backscatter color ratio for TTL Ci. The 
statistics are noisy because the sub-visible 
Ci are thin. Nonetheless, they provide 
information about this critical region when 
averaged appropriately. 
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Figures 3a and 3b: Cloud ice water path as 
a function of theta, and depolarization as a 
function of integrated backscatter, colored 
by 1064/532 nm color ratio (color scale). 

Figures 4a and 4b show the mean IAB and 
relative uncertainty for IAB and OD in very 
cold TTL clouds. Below -75 the uncertainty 
is dominated by the difficulty of detecting 
thin layers accurately. Note that we define 
“cold” clouds as having mid-layer 
temperature < -55°C, but when selected 
dynamically (theta > 355 K) the cloud mid- 
layer temperature < -65°C. 
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as a function of mid-cloud temperature, 
showing that clouds become thinner at cold 
temperatures in the TTL. Figure 4b shows 
the relative uncertainty for optical depth 
(red) and integrated backscatter (blue). 

4. CONVECTIVE INFLUENCE ON LOWER 
STRATOSPHERIC WATER VAPOR 

A model test, CALIOP space-based lidar 
and in situ aircraft observations all 
demonstrate that convective transport of 
water vapor and ice to the Tropical 
Tropopause Layer (TTL) and lowermost 
stratosphere is a significant process. 
Figures 5a and 5b show results from a 
simple cloud model that is added to a global 
domain-filling forward trajectory model 
(Schoeberl and Dessler, 2012) to test the 
impact of convectively lofted cloud ice and 
water vapor on the TTL and lowermost 
stratosphere. Results suggest that the anvil 
ice is needed to match CALIOP 
observations. 
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Figure 5: From Schoeberl et a!., 2014, 
shows the global distribution of CALIOP 
clouds for DJF 2008-9. The trajectory model 
is compared to the measurements, with and 
without MERRA anvil ice added to a simple 
cioud model. 



CALIOP ice water content for winter 2008-9 
shown in Figure 6 as an example of how 
much ice mass there is at high altitudes. 
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Figure 6: CALIOP zonally averaged ice 
water content during DJF 2008-9. Some of 
the apparent "noise" is due to less frequent 
overpasses at low latitudes. 


Measurements of cloud microphysical 
parameters in the TTL are few. Figure 7 
shows flight tracks from the CR-AVE 
mission out of Costa Rica in JF 2006. 


CR-AVE Field Mission 


CR-AVE = Costa Rican AURA 
Validation Experiment 


Platform: WB-57 
Instruments: 

CVI (just noise in SV Ci), Ozone, 
Chemistry, CVI, H20 vapor (!), 
Cloud/Aerosol Probes 

Cloud probe measurements: 

SPEC, Inc. - 2DS: 10 - 1400 urn 
DMT -CAPS: 1 - 1600 um 



Lawson analysis: ~10,000 sec of data 
This analysis: ~6,900 sec. of data 


2DS not in archive for: 
Transit flight 1, RF5, RF11 


Figure 7; Flight tracks from CR-A VE 



A SPEC, Inc. optical array probe (2D-S) 
measured the particle size distributions, 
extinctions and IWC in the TTL and 
lowermost stratosphere from the WB-57. 
IWC measurements are shown in Figure 8. 


CR-AVE Flights in 2006 
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Figure 8: SPEC 2D-S IWC composite 


Figure 9 shows the corresponding CALIOP 
cloud fraction climatology. 
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Figure 9: Cloud fraction, San Jose Costa 
Rica, DJF 2008-2012. 



An interesting result from CR-AVE is that 
there is more IWC at the cold point, but the 
largest particles occur above 380 K in the 
lowermost stratosphere. 
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Figure 10: CR-AVE 2D-S IWC vs. Extinction 
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Figure 11: Average particle size distribution 
in the TTL measured by the 2D-S 
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These size distributions agree well with the 
observations of de Reus et al., 2009, from 
Darwin during SCOUT. Typical CPI images 
of these particles are shown in Figure 12. 
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Figure 12: SPEC CPI images of ice particles 
in the lowermost stratosphere during CR- 
AVE. 


5. EVALUATING EXTINCTION 

COEFFICIENTS 

A new compilation of extinction coefficients 
measured using aircraft particle probes 
(details in Heymsfield et. al., 2014) has been 
evaluated, and statistics of the distribution 
are shown in Figure 14a, in blue. The mean 
extinction clearly decreases at colder 
temperatures, as will the associated ice 
water content. CALIOP data is shown in 
black for all of the data, and in red for 
constrained (QC=1) extinction solutions. 
Plot 14b shows that the constrained 
extinctions are not representative of the 
ensemble of CALIOP extinctions within each 
temperature range, and that at cold 


temperatures the CALIOP mean is larger 
than the mean for the probe data. 
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Figure 14a and b: In situ measured 
extinction coefficients (blue) compared with 
CALIOP retrieved extinctions - all (black) 
and only constrained solutions (red). 

Figure 15 shows that for temperatures 
colder than -55 deg C, constrained solutions 
are 20% or less of the data. 
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Figure 15: CALIOP distribution of extinction 
solution types as a function of temperature. 
Most cold Ci clouds are transparent to 
CALIOP and unconstrained (blue); these will 
be assigned an a priori value of the lidar 
ratio to solve for o. 


Figures 15 a and b show median integrated 
backscatter and extinctions as a function of 
temperature. The green line on plot d) is 
made by increasing the default lidar ratio 
from 25 to 32. 


Median IAB for topmost cloud layers, December 2009 



Temperature 


Nighttime CALIOP Median Extinctions, December 2009 

0.4 

0.3 

o 0.2 
o 

2 0.1 
0.0 

-80 -60 -40 -20 0 

Temperature 

Figure 15: a) IAB and b) Extinction as a 
function of temperature. 

6. COMPARISON WITH MR 

A matched set of CALIOP and MR data from 
December, 2009 (as in Parts II and II) is 
used for comparison between MR and 
CALIOP retrievals of cloud layer optical 
depth, effective diameter and ice water path. 
The layers were selected to optimize the MR 
retrievals with these characteristics: single 
(topmost) transparent (to CALIOP) cloud 
layers of MR scene types 21 , 30 and 24 (see 
Gamier 2013) with measured reference 
radiances. These layers were then further 
classified as: 1 - Cold: centroid 

temperatures < -55 deg C (30k samples); 

2 - TTL: cloud bases > 14.5 km (6k 
samples) 
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Figure 16a: Occurrence of layer geometric 
thickness as a function of temperature for 
"cold" ice clouds (T < -55 deg C) 
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Figure 16b: As for 16a, but TTL layers only. 

In the month of December, 2009 there were 
only 44 Ci layers with constrained CALIOP 
retrievals in the TTL. These data points are 
labeled with crosses on Figures 16b, 17b 
and 18b. CALIOP optical depth is in general 
lower than the equivalent MR optical depth 
because the unconstrained retrievals are 
assigned a lidar ratio of 25 that is too low for 
these thin clouds. But for layers with small 
integrated backscatter, the constrained 
optical depths appear to be much higher 
than the MR, suggesting that the retrieved 
lidar ratios are too high for weakly scattering 
cloud layers, or that the IIR optical depths 
are too low for these clouds. 
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Figures 17a and 17b: 2*1 IR 12.05 pm 

optical depth compared with the CALIOP 
optical depth (layer integrated extinction) for 
"cold" and "TTL" ice clouds, respectively. 


We evaluated effective diameter as it is 
defined in Mitchell (2010): 

Deff=(3/ P ice)(IWC/a) (2) 

In this case we are using integrated 
quantities for the cloud layers, so the ratio is 
ice water path (IWP) to optical depth. 
CALIOP effective diameters are larger than 
those from the IIR, which may indicate that 
the Version 3 IWC parameterization is too 
large for colder temperatures. New analysis 
in Heymsfield et al., 2014 suggests updated 
IWC parameterizations that are currently 
being evaluated for possible use in CALIOP 
Version 4. 
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Figures 18a and 18b: IIR vs. CALIOP cloud 
ice water path for "cold" and "TTL" ice 
clouds. 

Figures 18 a and b show that the IIR and 
CALIOP ice water path for cold and TTL Ci 
clouds agrees fairly well, with CALIOP 
slightly higher even while having low optical 
depths, so the differences apparently cancel 
each other somewhat. The CALIPSO team 
is currently evaluating lidar ratios and iwc 
parameterizations for the new Version 4 
Level 2 product. 


The CALIPSO CALIOP and IIR Version 3 
products are available at: 

NASA LaRC: ( http://eosweb.larc.nasa.gov/ ) 

I CARE: (http://www.icare.univ-lille1.fr/ ) 
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